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Abstract - This paper presents a unified methodology of combining heuristic fuzzy design and FEM verification to design 
optimally vessels with a branch joint. This method is applicable to steels and composites. Background for this study is the 
recognized need for constructing vessels with branch tubes needed for processing liquids and gases with minimal ecological 
and maintenance problems. The methodology of fuzzy optimum design is used. The goals and constrains are expressed in a 
consistent formulation. First, design variables like wall thickness are defined discretely within ranges. Then decision 
variables are formulated like cost and safety factors. The total goal is maximization of the end user satisfaction on the 
design. It is product of satisfaction functions of decision variables. The stresses are calculated by reasonable free body 
models and notch factors. Two steels are studied, a basic low strength steel (LS) and a high strength steel (HS). At low 
pressure p=0.1 MPa the LS vessel is 4.2 times more satisfactory than the HS vessel. At higher pressure p=l the LS vessel is 
0.4 times less satisfactory than HS vessel. An analytical stress result agrees reasonably with FEM results at a test pressure. 
The optimal choice depends not only on economics and technology but also on the societal and environmental changes and 
megatrends. This methodology can be used to explore novel concepts. 

Keywords: pressure vessel, fuzzy optimum design, steels, fatigue. 

I. Introduction 

This study is motivated by the recognized global need for constructing vessels with branch joints needed for processing 
liquids and gases with minimal ecological problems. Cylindrical vessels with elliptical ends are reasonable choices. 

Customers require a reliable lifetime of tens of years with minimal initial cost, maintenance and still stand costs. These 
goals can be formulated and solved by applying heuristic fuzzy optimum design with FEM verification. 

There are many choices for vessel materials. All materials are composite from macro to microstructural level. In certain 
corrosive loadings fiber reinforced composites (FRP) are a competitive choice to highly alloyed steel and in others not. In 
this task a reasonable selection of steel and composite materials are options. 

Composites for structural engineering are commonly made by combining materials from five commonest material groups, 
metals, ceramics glasses, plastics and polymers. Even only two components one can obtain many useful combinations of the 
properties. The one with largest volume fraction is called matrix. Others can be in form of fibres, platelets and globular 
forms. 

The fibre reinforced plastics are designed to combine the strength of strong fibres in desired direction, chemical, static and 
fatigue strength, impact, environmental and thermal endurance. Composite material design is discussed by [1] (Agarwal, 
Broutman, 1990), by [2] (Barbero, 1990) and [3] (Swanson, 1997). A novel generalized failure criterion is proposed by [4] 
(Knops, 2008). Now a modified criterion is developed. 

Basic vessel theory and designs of steel materials are discussed by [5] (Harvey, 1991). The main method of vessel design is 
the basic shell theory. Mandatory design rule codes are based on theory, tests. FEM is used to verify the analytical results. 

The vessel design goal can be expressed easily as fuzzy satisfaction of the end user. This method is based on results of [6] 
(Diaz, 1988). [7] (Martikka, Pollanen, 2009) have applied multi -objective optimisation by technical laws and heuristics to 
vessel design. This same methodology is applicable to a wide range of optimization tasks. Design of metal and non-metal 
materials for pressure vessels is discussed by [8] (Martikka, Taitokari, 2011). 
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In these tasks the reliability based design can be applied as discussed by [9] (Leitch R.D, 1988), and [10] (Dhillon, and 
Singh, 1981). Optimum design of FRP vessels is considered by [11] (Martikka, Taitokari 2006) under variable pressure 
loads. Design of FRP vessels with internal dynamical load due to mixing rotor is discussed by [12] (Taitokari, Martikka, 
2009). External dynamic load to vessels are often seismic and wind loadings. Optimal design of seismically loaded vessels is 
discussed by [13] (Martikka, Pollanen, 2011). Endurance dimensioning of steels are discussed by [14] (Gurney, 1978) and 
[15] (Meyer, 1985). 

Design of welded tube to vessel joint for liquid gas is discussed by [16] (Martikka, Taitokari, 2012). Welded beam joint 
design is discussed by [17] (Martikka, Pollanen, 2008). Optimal dynamical and statics design of steel bridges is discussed by 
[18] (Martikka, Taitokari, 2013). 

The main purpose of this study is to develop methodology by which branched steel vessels can be optimised using fuzzy 
optimum design with goal of maximization of user satisfaction. 

II. Design approach 

The objective in this section is to present an overview of the present approach for obtaining a satisfactory design. Composites 
offer more design variables than steels allowing a trade-off use of many decision variables, like low cost, low weight, high 
safety at extreme loadings. 

The goalsof designing of thin shelled structures for process industry 

Successful producers are able to recognize the needs of customer . The top level goal is psychological, to obtain maximal 
satisfaction of the customer on the total costs and reliable time of utilization of a product. The task of the designer is to 
devise a method and apply it to reach the goal. 

A. Satisfaction on technology of the design concept 

This is mainly decided by the design engineers. Technical decision variables are factors of safety. 

B. Satisfaction on cost is mainly defined by the customer. 

It includes the total cost, maintenance and recycling and ecology costs. The total goal is maximal satisfaction on technology 
and total life cycle costs. 


III. The OBJECT OF DESIGN 

The objective in this section is to present an overview of the object of the present case study. 

3.1 A typical vessel structure 

In Figure 1 typical structures and loadings are shown. 



Figure 1: An overview sketch of vessel; a) dominant stresses and main radii; 
b) Bottom designs; c) Surrogate modeling of internal stresses 


The basic geometry of the structure is shown in Figure 2a. The typical features are support, bottom opening, and top opening 
and side branch opening. The main load is now internal pressure p inside. Branch tube may be closed by two options: a 
conventional model and a quick closure model. The critical locations are shown. 

Reasonable dimensions are:Geometry of the main vessel A : R = 1 m, 77 = lm, tube branch r = .5m, / = lm. Reasonable 

choices for the p = IMPa case are T = 0.013m, f h = 0.005m, 7] = 0.013m 777) = 1.2,777,, = 1 whence 7]= 771.2 = 0.013/1.2, 
L r = T= 0.013m. 
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A typical industrial vessel is shown in Figure 2b. 




<b) 


Figure 2: a) The vessel in this case study. 

B) A TYPICAL INDUSTRIAL VESSEL AND FEM MODEL IN THIS STUDY. 


3.2 Joint stresses 

The complex bending stresses at the joint can be estimated using a simple surrogate bent beam model shown inFigure 3 
schematically. The corresponding FEM model is shown in Figure 4. 



Figure 3: A simple bent beam model for estimating the joint stresses . 


The bending force on the model is F b , It is one quarter of the resultant pressure force acting on the branch head plate. It 
causes a bending moment. The lever is some radial flange distance A r. The width of the beam is peripheral length around 
the tube B of thickness f h . 
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Substitution gives the nominal bending stress of B at joint 
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Here x = extent of stiffening influence by this fictive model. 

The axial nominal stress at the shell B is caused by the pressure 
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1 MPa ■ 0.5 m 
2 • 0.010m 


= 25 MPa 


The total summed nominal axial stress in tube B is sum 


S x — S x „ + S xb — ct Bx 


f . \ 

1 + 6 — 

t h j 


= C 7 , 
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The nominal stress in the y direction of the branch B is 
Sy ~ & Bt ~ 2 <T Bx 
The equivalent Mises stress at B is 


ct VA b = yjs x 2 + s y 2 -s x s y = o B J( l + 6x) 2 +4-2-(l + 6x) = o Bx G(x) 

The peak stress is obtained by multiplying it by a stress concentration factor K , 

^VAB max ~ ^t^VAB 

This is a heuristic tentative model. The following rough estimates are feasible 


(3) 

(4) 

(5) 

(6) 

(7) 
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The G is a rough model only 

x — 0.2 => G « 2 , x « 0.4 => G « 3 , x « 0.6 => G « 4 

The stress concentration is estimated in the program using choices K t = 3 and G =3 


^YABmax = K x<*vxb = ' G ' ^x B = ^totai ■ ^xB = 3 ‘ 3 • 25 MPa = 225 MPa 


( 8 ) 


(9) 


The stress distribution with FEM is shown in Figure 9b. The notch is sharp giving a high stress concentration of about 10 
since this is a linear elastic FEM element model. 



Figure 4: The Mises stress by FEM [20] at the joint. 

FEM gives maximum value of Mises stress as 1045 MPa. 

The stress concentration by FEM is ATfem =10 due to a very small notch radius 

a \w m Jy FEM ’ shar P) = K ¥EU -cr xB * 10*100 MPa = 1000 MPa (10) 

If the FEM model notch geometry is changed to have a moderate notch radius then the stress concentration will be 
acceptably small 


°"vABmax ifEM > smooth ) = K FEM • cr xB ~ 3 • 1 00 MPa = 300 MPa 


( 11 ) 
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In the analytical model the total stress concentration is a heuristic estimate, 

(7 VA t> = K . ■ G • o- xR = 3 ■ 3 • 25 MPa = 225 MPa 4 ■ 4 • 25 MPa = 400 MPa 

VAb max t xb 

The FEM and analytical stress peaks are fairly close. 

The mean Mises stress and amplitude Mises stresses are 

^Vm = ^VAB.max ’ ^Va = 0.2o" VlT1 , AcJ = 2<Ty a 


( 12 ) 


(13) 


Here the assumption is made that the amplitude stress is 0.2 times the mean stresses. These are used in fatigue life estimation 

3.3 Simplified analysis of stresses in a tube branch connection 

The objective in this section is to present an overview of the present approach to get the tube branch stresses using a 
simplified analysis based on FBD balance shown in Figure 5. 



Figure 5: Simplified analysis of stresses in a tube branch connection; 
a) sketch; b) Notch factor model 


The pressure force is balanced by the resultant force due to shell tangential stress. 
The nominal tangential stress acting on the cross sectional area of the nozzle tube is 

'A„ 


<?B=P 


B 


A=^HR + lr,B =ll + (\H- r)r, <7 B , m „ = K,a„ 


(14) 


Here the stress areas are A p and/7. c Bma x is the maximal stress and K t is the notch factor . 

Here the notch factor is estimated with an approximation of the model in [19] (Dubbel, 1994). 
Flat plate surrogate model is used to estimate the stress concentration factor. 

Two options are given for model parameter q , <:/=(). 5 for bending q= 1 for tension. 


K 


t, stepped, plate 




The tube wall may be about f h = 5. .8 mm. Then a= 0.5 f h . Radius may be p = 1mm. The estimated range is 

K, = l + r/. 


4/77/77 


t,. stepped plate 


1/77/77 


= 1 + 2q = 2(q = 0.5). ...3(4 = 1) 


(15) 


(16) 


IV. Fuzzy goal formulation using decision variables 

The objective in this section is to present of an overview of the present approach to get a satisfactory design. The design 
variable vector x = (load functions, geometry, materials) is the tool. 

Using them the decision variable event vector s = (cost, factors of safety...) = s(x) is formulated. They are defined fuzzily 
(Diaz, 1988). The total event is decision variable s and it is intersection of the chosen decision variables ,y k 
s = ,V| 


( 17 ) 
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The design goal is maximization of the total satisfaction of the customer on the product 

P(s) => p(s) = P(s , ) • P{s 2 ) • • P{s n ),Q = max P 

Now all goals and constraints are formulated consistently by one flexible fuzzy function. This is illustrated in Figure 1 1 


(18) 


V. Design variables 

The aim in this section is to present the main steel options with appropriate property data and estimations. Composite 
materials can also be activated with their appropriate criteria. 

5.1 Material design variables 

Steel option design variables are shown in Table 1. The basic steel S52 is commonly used in steel vessel constructions. The 
OX steel is used when high yield strength is needed. The stainless steels 904L is used when corrosion resistance is needed. 
An alternative to it is FRP. The E25Cr is used in corrosion resistant castings. 

Table 1. 

Steel material variables. The Paris law C and m parameters are calculated with Gurney’s 


model[14 

(Gurney, 1978). Stress rations = ctmin/ cj M in = o, K m (N-mm 3/2 ), 


OXsteel im =6 

St52 steel im=7 

Stainless 904L im=3 

E25Cr steel im=4 

Yield strength (MPa) 

Re(6) = 1000 

Re(7) = 335 

Re(3) = 220 

Re(4) = 500 

Tensile strength (MPa) 

Rm(6) = 1100 

Rm(7) = 500 

Rm(3) = 490 

Rm(4) = 700 

unit cost (eur/kg) 

Cm(6) = 20 

Cm(7) = 5 

Cm(3) = 5 

Cm(4) = 5 

density (kg/m 3 ) 

rho(6) = 8000 

rho(7) = 8000 

rho(3) = 8000 

rho(4) = 8000 

ecological value 

eco(6) = .5 

eco(7) = .5 

eco(3) = .7 

eco(4) = . 1 

corrosion resistance 

corres(6) = .5 

corres(7) = .5 

corres(3) = .8 

corres(4) = .9 

Elastic modules (MPa) 

E(6) = 205000 

E(7) = 205000 

E(3) = 205000 

E(4) = 205000 

Threshold intensity(Nmm" 3/2 ) 

Kth(6) = 275 

Kth(7) = 190-144 -Rs 

Kth(3) = 3MPam‘ /2 

Kth(4) = 3MPam‘ /2 

Initial crack size (mm) 

a0(6) = 1 

a0(7) = 1 

a0(3) = 1 

a0(4) = 1 

Paris C parameter 

C(6) = 4.64E-12 

C(7) = 1.67E-14 

C(3) = 1.67E-14 

C(4)= 1.67E-14 

Paris m exponent 

m(6) = 2.52 

m(7) = 3.36 

m(3) = 3.36 

m(4) = 3.36 


5.2 Functional design variables and parameters 

The pressure p was now one input parameter. It can be changed to a variable if optimal loading is sought. 

5.3 Geometrical design variables 

Here the design variables related to the vessel and branches are 

T= wall thickness with discrete values in the range 0.005. . .0.013m, wall thickness h of the branch tube with range 
0.005. ..0.013 and relative height c/R with range 0.45. ..0.65. The end wall thickness ratio T/T wassettol. 


T = T(iT),iT = 1...8 h = f/?(ith),ith = 1...8, xci?(ixcR) = — ,ixcR = 1...8 

R 


The design variables related to the stiffener rings are. E r = E, for FRP and E r = E for steels 
The ranges are for x^O.5.,.1.4 and forx 2 =1...2 and for x 3 = 1..2.5 


x l - — = A'£>(ixEr), 
E. 



T T 
L T 

r r 


xAr( ixAr),x 3 



(19) 


( 20 ) 


VI. Decision variables for defining goals 

The objective in this section is to present of overview of the decision variables used. 
Decision variables s(k) depend on the design variables x(i) , s = s(x). 
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The results for the optimal low alloy vessel at load p =1 MPa are the following: The optimal geometry is T =0.013m, t\, = 
0.01, r e = 0.013, xCR = c/R = 0.45, xeT = 1, xEr = 1, xAr = 2 . Satisfactions P(s k ) and s k are obtained for the optimal low 
alloy construction 

6.1 Goal of obtaining low cost K or decision variable s. 

Cost is now the material cost of the shells and plates. 

6.1.1 The ellipsoidal end of the vessel 

The bottom radius is the same as the radius of the main tube A. Now R is given and height c of the ellipsoid is a design 
variable. The radii of curvature of the shell are 

2 

c 

r 2 =ci = b = R , r x = — , c = xcRR ( 21 ) 

R 

Radius of curvature at the calotte at x = 0, z = c 



Figure 6: Geometry of the ellipsoidal calotte end of the vessel 


The total cost the first decision variable. Same unit cost C is assumed for all components 

5, = K = C(m A + 2 m e + m B + m Be ) 


(23) 


6.2 Goal of obtaining satisfactory safety factor for the main shell in tangential direction 


Tangential direction is chosen as the strong direction L. The branch opening raises the stress 


D 

V - N - _ml 

Jv tA ’ 


<T, 


pRH pR 1 


a, 


tA 


tA 


(H~r)T 


— > 


T 1 - 


= cr, 


H 


6.3 Goal of obtaining satisfactory safety factor for the main shell A in axial direction 


(24) 


The decision variable is now the safety factoragainst axial stress o xA 


= ^xA = 


/? 


xA 


= 


p\jiR 2 


TR(±7t-<f>) 4 T 


= \^- = < t t 4 = f (25) 


6.4 Goal of obtaining satisfactory safety factor for the tube shell B in tangential direction 

Tangential direction is now chosen as the strong direction Land the load stress o tB are 
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s 5 ~ — 

6.6 Goal of obtaining satisfactory safety factor for the calotte of the main shell 




, <7 lB (J [5 


(B 


pr 

t 


— (26) 


Safety factor of the calotte of the calotte of the main shell A 

R «* r 

’ ' ( 


s — yv 

13 7 iV calot 


G 


calot 


T -T 

’ calot * e ’ *7 


calot 


calot 


P r a„ot 

T, 


(27) 


calot 


Here R m is UTS of FRP in x direction, in steel it is f? m (im) 


6.6 Goal of obtaining satisfactory safety factor for joint of the shells based on Mises equivalent stress 

The decision variable is the safety factor 

Rjjm ) 

JV AB 

Cr VAB,max 

6.7 Goal of obtaining satisfactory safety factor for joint of the shells using the simple method 

The decision variable is 

n 

s = N = 23 

Jv xB 


xB,max 


The maximal stress is 


^xB.max ^t^B’^t * ^ 


6.8 Goal of obtaining satisfactory safety factor for the main shell maximal stress location 

The decision variable is factor of safety 


*10 — ^pert.A 


R 


or 


V , pert A 

The nominl stress S and the perturbation stress acting on vessel A are 

1-1 




(28) 


(29) 


(30) 


(31a) 


2 T 2-0.013 


= 38 MPa, r l =0.3R , e = 0.146, cr VpeitA = 1.15 (31b) 


6.9 Goal of obtaining satisfactory safety factor against buckling of the main shell 


The pressure loading at low pressure causes a buckling risk. 

To avoid buckling risk the tangential stress should be tensile and the factor of safety within a range 


S 


li 


= N 


buckl,A 


=0.605 E- 
a 2 R 


Her scr is The buckling strength for a compressed ideal cylinder with wall T and radius R 
The principal stresses are 


03 = < 7 , = 


pR 

zP 


® 8 ^2 


pR 


1- 


R 

2c 2 


2 A 


pR 

T. 


R 

2r, 


i y 


(32) 


(33) 


6.10 Goal of obtaining long enough crack initiation life at the joint 

The decision variable is 

*,2=^-. 2 < ^ l2 < 20000 


(34) 
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Satisfaction on s 12 is P(s 12 =6.8 ) = l.Now it is by passed by setting P = 1 but can be activated. 

This method of calculating fatigue life A life combines the Haigh diagram of modified Goodman type and the S-N diagram 
according to (Meyer, 1985). Its is graphically shown in Figure 7 


N life = 10 a ,A = log 


v t v e 

3 

_a-vjc 2 _ 

_log(V e / c)_ 


w ^"va w 'G m 

’ a ” R. ’ m ~~ 


,V e = — ,c = 0.9 (35) 

R 


HereV a is relative effective stress amplitude, V m is relative effective mean stress and V e is relative effective corrected fatigue 
strength. Equivalent stress at the joint of A B at the upper point of the joint. The mean Mises stress and amplitude Mises 
stress are 


<J 


V.AB 


cr, " + a~ + 3 1 


LT 


1 

• 2 


^Vm = <?■ 


V.AB 


fT Va = 0.2cr Vm 


(36) 



Figure 7: The method of calculating fatigue lives of crack initiation time from normal 

MEAN STRESS AND AMPLITUDE STRESS VS. S-N DIAGRAM 


6.11 Goal of obtaining satisfactory fatigue crack propagation life 

The decision variable is now the Paris crack propagation life prediction 


N 


S 13 


pans 


10 6 


■O ' 1 = Gamin < G 3 < Gamax = 100 °Q 


(37) 


The fatigue life in number of cycles from initial to final crack length is 


N(im ) = 


1 


C(imj\ /»(//»)- 1 )(act v 


YyJX 


i(im) 


:/„(//»)- 


7-1 


: (lm) 2 


1-1 


(38) 


6.12 Goal of obtaining satisfactory crack stress threshold safety factor 

The decision variable is now safety factor against crack stress threshold exceedance, ao is initial crack size and a th is the 
threshold crack size .Now satisfaction was set to 1 


s = N = 

J 14 JY th 


«th 

mm 


a o ^ 

[imj 

mm] 


1 

n 


A K 
YAa 


th 


max J 


2-[MPam‘ /2 

2-100[MPa 


0.03mm 


(39) 


6.13 Goals of obtaining satisfactory ecology and corrosion resistance 

Goal of obtaining satisfactory ecological and corrosion resistance value depends on the material selection code im 

5 16 = eco(im), 0.05 < s 16 < 0.95, ^ = corredjm), 0.05 < ^ < 0.95 (40) 


6.14 Goal of obtaining satisfactory factor of safety for the ring stiffener 


Page | 14 
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The decision variable is factor of safety against stiffener ring yielding depends on the UTS of the material 

s K =N mi =-^, R,„=R„ (im) 

^"max, nng 

A model of a ring stiffener on a cylinder under internal pressure is shown in Figure 8. 



(41) 


Figure 8: Ring stiffener on cylinder stiffener model. 


VI. Results and discussion 


The aim in this section is to present the main results are depending on the pressure and material selection. The total 
satisfaction P is product of two satisfactions 


P - P P P 

tech cos t tech 



(42) 


Figure 9 shows comparison of the results for the low and high strength steel vessels depending on pressure p. 
Here R = Ratio of technical satisfactions ,R= 0.1(P tech (LS) /P, ec h(HS) ) 



1 

P cost 


Cost/10 


0.5 


pressure 


A) 



1 

P cost 

Cost/20 

0.5 


0 


B) 


pressure 
P ,( MPa) 


FIGURE 9: TOTAL SATISFACTION AND COST SATISFACTION AND COST OF THE VESSEL 
DEPENDING ON THE PRESSURE P. A)LOW STRENGTH LS B) HIGH STRENGTH HS 


VII. Use of FEM to verify optimal concepts 


The results of FEM are shown in Appendix Al. 
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The strategy of using FEM in design was as follows. First the optimal concept to be tested was created, optimized and 
selected 

Then the stresses and strains and buckling modes are calculated. 


VIII. Conclusions 

The main conclusions of the present study can be summarized as follows 

A§ The motivation to this study is to develop a methodology for obtain optimal concepts vessels with tube 
branches 

B§ Thedesign variables are 

• The material options are metal and non-metals. 

• The geometrical options are wall thicknesses 

C§ Using the methodology of unified fuzzy design the goals and constraints are defined by the same formulation. 

• Satisfactory cost and factors of safety and other decision variables are obtained using few discrete design variable 
selection options 

D§ Next satisfaction functions are defined for each decision variables. Then the design goal is to maximize their 
product. 

E§In the present study the low strength and high strength steels are competitive choices with strong and weak properties. 

The total satisfaction was product of costs satisfaction and technical satisfaction. 

The total satifaction for LS steel is P tot (LS)= P cost (LS) P tech (LS) 

The total satifaction for HS steel is P tot (HS)= P cost (HS) P tech (HS) 

At low pressure p=0. 1 MPa P tot (LS)/ P tot (HS) = 4.2: 1 
At high pressure p=l MPa P tot (LS)/ P too (HS) = 0.4 :1 

Technical satisfaction behaviours differ with HS and LS vessels 

• At low pressure p = 0. 1 the ratio P tech (LS)/ P te ch(HS) = 2.69/0.86 = 3.1:1 

• At moderate p = 0.5 and at p=l the ratio P tech (LS)/ P, ec h(HS) = 0.5/1 .7 =0.3:1 

A possible explanation is the following 

• At low pressure p = 01 MPa the strength properties of HS are underutilised and satisfaction is low 

• At higher pressures the good strength properties of HS can be better utilized by using thicker walls 

F§ the optimality of the concept obtained by analytical optimization is verified by FEM models. 

G§The present methodology will be used in future to explore new innovation concepts of various materials which are needed 
in the near future. 
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APPENDIX Al. FEM RESULTS 

A FEM model was made to verify the optimality of the low strength steel vessel at p = IMPa inner pressure. The FEM 
element model is shown in Figure Al-1. 



a) b) 

FIGURE Al-1. A) FEM MODEL SHOWING THE ELEMENTS. THERE ARE 63739 ELEMENTS 
AND 63624 NODES [20]. B) THE MAXIMAL STRESSES ARE AT THE UPPER AND LOWER BRANCH JOINTS 

The geometry is nearly the same as in the analytical optimal model. 

• vessel wall T= 0.013, 

• branch tube wall r h = 0.010 in analytical model , f h = 0.0050 in FEM model 

• calotte wall T c = 0.013, 

In the analytical model 

• vessel radius R= lm half height H=lm, 

• branch tube radius r = 0.5m , length 1 = lm 

• vessel calotte relative height xcR =c/R =0.45, calotte relative wall xET= TJT= 1 , 

• stiffener elastic modulus xEr = E/E r =l, stiffener area x A R = T 2 /A r = 2 
In the FEM model nearly the same ends are used 

Lower an upper elliptical ends are according the Korbbogen code. The peripheral stress in the branch tube is 103 by FEM 
and 100 by analytical model. 


